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ABSTRACT 



^Observations reported here were obtained at the MMT Observatory, a joint facility of the University of Arizona 
and the Smithsonian Institution. 



- 2 - 



We have carried out an L' and M band Adaptive Optics (AO) extrasolar planet 
imaging survey of 54 nearby, sunlike stars using the Clio camera at the MMT. Our 
survey concentrates more strongly than all others to date on very nearby F, G, and K 
stars, in that we have prioritized proximity higher than youth. Our survey is also the 
first to include extensive observations in the M band, which supplemented the primary 
L' observations. These longer wavelength bands are most useful for very nearby systems 
in which low temperature planets with red IR colors (i.e. H — L' , H — M) could be 
detected. The survey detected no planets, but set interesting limits on planets and 
brown dwarfs in the star systems we investigated. We have interpreted our null result 
by means of extensive Monte Carlo simulations, and constrained the distributions of 
extrasolar planets in mass M and semimajor axis a. If planets are distributed according 
to a power law with dN oc M'^a^dMda, normalized to be consistent with radial velocity 
statistics, we find that a distribution with a = —1.1 and f3 = —0.46, truncated at 110 
AU, is ruled out at the 90% confidence level. These particular values of a and /3 are 
significant because they represent the most planet-rich case consistent with current 
statistics from radial velocity observations. With 90% confidence no more than 8.1% of 
stars like those in our survey have systems with three widely spaced, massive planets 
like the A-star HR 8799. Our observations show that giant planets in long-period orbits 
around sun-like stars are rare, confirming the results of shorter- wavelength surveys, and 
increasing the robustness of the conclusion. 

Subject headings: planetary systems, techniques: IR imaging, intrumentation: adaptive 
optics, astrometry, binary stars 



1. Introduction 



Nearly 400 extrasolar planets have now been discovered using the radial velocity (RV) method. 
RV survey s currently have good s tatistical coinpleteness only for planets with periods of less than 
ten years (jCumming et al.l l2008l : iButler et al.l l2006l : iFischer Sz Valentil l2005l ) , due to the limited 
temporal baseline of the observations, and the need to observe for a complete orbital period to 
confirm the properties o f a planet with con fidence. The masses of discovered planets range from 
just a few Earth masses (jBouchy et al.ll2009l ) up to around 20 Jupiter masses (Mjup). We note that 
a 20 Mjup object would be considered by many to be a brown dwarf rather than a planet, but that 
there is no broad consensu s on how to define t he upper mass limit for planets. For a good overview 



of RV planets to date, see iButler et al.l (|2006l ) or http://exoplanet.eu/catalog-RV.php 



The large number of RV planets makes it possible to examine the statistics of extrasolar planet 
populations. Several groups have fit approximate power law distributions in mass and semimajor 
axis to the set of known extrasolar planets (see for example ICumming et al.l (120081 )). Necessarily, 
however, these power laws are not subject to observational constraints at orbital periods longer 
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than 10 years - and it is at these orbital periods that we find giant planets in our own solar system. 
We cannot obtain a good understanding of planets in general without information on long period 
extrasolar planets. Nor can we see how our own solar system fits into the big picture of planet 
formation in the galaxy without a good census of planets in Jupiter- and Saturn-like orbits around 
other stars. 

Repeatable detections of extrasolar planets (a s opposed to one-time mic rolensing detections) 
have so far been made by transit detection fe .g. ICharbonneau et al.l ([20001)), by RV v ariations 
dMavor k Quelojll995l ). by astrometric wobble (jBenedict et al.ll2006l ) , or by direct imaging ([Marois et al, 
20081 ). Of these methods, transits are efficient only for detecting close-in planets. As noted above, 
precision RV observations have not been going on long enough to detect more than a few planets 
with periods longer than ten years, but even as RV temporal baselines increase, long period planets 
will remain harder to detect due to their slow orbital velocities. The amplitude of a star's astro- 
metric wobble increases with the radius of its planet's orbit, but decades- long observing programs 
are still needed to find long-period planets. Direct imaging is the only method that allows us to 
characterize long-period extrasolar planets on a timescale of months rather than years or decades. 

Direct imaging of extrasolar planets is technologically possible at present only in the infrared, 
based on the planets' own thermal luminosity, not on reflected starlight. The enabling technology 
is adaptive optics (AO), which allows 6-lOm ground-based telescopes to obtain diffraction lim- 
ited IR images several times sharper than those from HST, despite Earth's turbulent atmosphere. 
Theoretical models of giant planets indicate that such telescopes should be capable of detecting 
self-luminous giant planets in large orbits around young, nearby stars. The stars should be young 
because the glow of giant planets comes from gravitational potential energy converted to heat in 
their formation and subsequent contraction: lacking any internal fusion, they cool and become 
fainter as they age. 

Several groups have published the results of AO imagin g surveys for extrasolar p lanets around 
F, G, K, or M star s in th e last five years (see for example 



Masciadri et al 



diooi); iBiller et al.l ^200^ ): iLafreniere eTZI (|20o3); and Ichauvin et al 



j2005|); iKasper et al 
(l20inl )). Of these, most 



have used wavelengths in the 1.5-2.2 /xm r ange, corresponding to the astronomical H a nd Ks filters 
(jlVIasciadri et al.l boosl : biller et al.] boO?! : ILafreniere et al.l booj : Ichauvin et aD boiol ) . They have 
targeted mainly very young stars. Because young stars are rare, the median distance to stars in 
each of these surveys has been more than 20 pc. 

In contrast to those above, our survey concentrates on very nearby F, G, and K stars, with 
proximity prioritized more than youth in the sample selection. The median distance to our survey 
targets is only 11.2 pc. Ours is also the first survey to include extensive obser vations in the M ban d, 
and only the second to search solar- type stars in the L' band (the first was iKasper et al.l (j2007l )). 
The distinctive focus on older, very nearby stars for a survey using longer wavelengths is natural: 
longer wavelengths are optimal for lower temperature planets which are most likely to be found in 
older systems, but which would be undetectable around all but the nearest stars. More information 
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on our samp l e selec tion, observations, and data analysis can be found in our Observations paper, 
Heinze et alj (|20ld ). which also details our careful evaluation of our survey's sensitivity, including 
extensive tests in which fake planets were randomly placed in the raw data and then recovered 
by an experimenter who knew neither their positions nor their number. Such tests are essential 
for establishing the true relationship between source significance (i.e. 5a, 10a, etc.) and survey 
completeness. 

Our survey places constraints on a more mature population of planets than those that have 
focused on very young stars, and confirms that a paucity of giant planets at large separations from 
sun-like stars is robustly observed at a wide range of wavelengths. 

In Section [21 we review power law fits to the distribution of known RV planets, including the 
normalization of the power laws. In Section [3l we present the constraints our survey places on the 
distribution of extrasolar giant planets, based on extensive Monte Carlo simulations. In Section [J] 
we discuss the promising future of planet-search observations in the L' and especially the M band, 
and in Section [5] we conclude. 



Statistical Distributions from RV Planets 



Nearly 400 RV planets are known. See iButler et al. for a useful, conservative list- 

ing of confirmed extrasolar planets as of 2006, or http : / / exoplanet . eu/ catalog- RV . php for a 
frequently-updated catalog of all confirmed and many suspected extrasolar planet discoveries. 

The number of RV planets is sufficient for meaningful statistical analysis of how extrasolar 
planets are distributed in terms of their masses and orbital semimajor axes. The lowest mass 
planets and those with the longest orbital periods are generally rejected from such analyses to 
reduce bias from completeness effects, but there remains a considerable range (2-2000 days in 
period, or roughly 0.03-3.1 AU in semimajor axis for solar-type sta rs; and 0.3-20 Mjup in mass) 



where RV searches have good completeness (jCumming et al 



20081 ). There is evidence that the 



shortest period planets, or 'hot Jupiters,' represent a separate population, a 'pileup' of planets in 
very cl ose-in orbits that does not follow the same statistical distribution as planets in more distant 
orbits (jCumming et al.l 120081 ). The hot Jupiters are therefore often excluded from statistical fits 
to the overall populations of extrasolar planets, or at least from the fits to the semimajor axis 
distribution. 



Gumming et al.l (120081 ) characterize the distribution of RV planets detected in the Keck Planet 



Search with an equation of the form 



dN = CoM"^P^^dln(M)dln(P). 



(1) 



where M is the mass of the planet, P is the orbital period, and Cq is a normalization constant. 
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They state that 10.5% of solar-type stars have a planet with mass between 0.3 and 10 Mjup and 
period between 2 and 2000 days, which information can be used to derive a value for Cq given 
values for the power law exponents ol and They find that the best-fit values for these are 
ai = —0.31 ± 0.2 and /3l = 0.26 ± 0.1, where the l subscript is our notation to make clear that 
these are the exponents for the form using logarithmic differentials. 

In common with a number of other groups, we choose to represent the power law with ordinary 
differentials, and to give it in terms of orbital semimajor axis a rather than orbital period P: 



dN = CoW^a^dUda. (2) 

Where Co, of course, will not generally have the same value for Equations [J and [2j Manipu- 
lating the two equations and using Kepler's Third Law makes it clear that 



a 



OiL 



(3) 



and 



/3 = ;^/3i 



1. 



(4) 



The ICumming et al.l (j2008l ) exponents produce a = —1.31 ± 0.2 and /? = —0.61 it 0.15 when 
translated into our form. The mass power law is well behaved, but the integral of the semimajor axis 
power law does not converge as o — > oo, so an outer truncation radius is an important parameter 
of the semimajor axis distribution. 



Butler et all (j2006l ) presents the 2006 Catalog of Nearby Exoplanets, a carefully described 
heteroge nous sample of plan ets detected by several different RV search programs. With appropriate 
caution, I Butler et al.l (|2006l ) refrain from quoting confident power law slopes based on the combined 
discoveries of many different s urvey s with different detection limits and completeness biases (in 
contrast, the ICumming et al.l (|2008l ) analysis was restricted to stars in the Keck Pl anet Search 



which were uniformly observed up to a given minimum baseline and velocity precision). iButler et al 



(|2006l ) do tentatively adopt a power law with the form of Equation [2] for mass only, and state that a 
appears to be about -1.1 (or -1.16, to give the exact result of a formal fit to their list of exoplanets). 
However they caution that due to their heterogeneous list of planets discovered by different surveys, 
this power law should be taken more as a descriptor of the known planets than of the underlying 
distribution. They do not quote a value for the semimajor axis power law slope /3. 



Based mostly on Gumming et al. ( 2008 ). but considering Butler et al. (j2006l ) as helpful addi- 
tional input, we conclude that the true value of the mass power law slope a is probably between -1.1 
and -1.51, with -1.31 as a good working model. The value of the semimajor axis power law slope 
/3 is probably between -0.46 and -0.76, with -0.61 as a current best guess. The outer truncation 
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radius of the semimajor axis distribution cannot be constrained by the RV results: surveys hke 
ours exist, in part, to constrain this interesting number. 

The only other result we n eed from the RV searc hes is a normalization that will allow us to 
find Cq. We elect not to use the ICumming et alJ (j2008l ) value (10.5% of stars having a planet with 
mass between 0.3 and 10 Mjup and period between 2 and 2000 days), because this range includes 
the hot Jupiters, a separate population. 



We take our normalization instead from the Carnegie Planet Sample, as described in lFischer &: Valenti 



(|2005l ). Their Table 1 (online only) lists 850 stars that have been thoroughly investigated with RV. 
They state that all planets with mass at least 1 Mjup and orbital period less than 4 years have 
been detected around these stars. Forty-seven of these stars are marked in Table 1 as having RV 
planets. Table 2 from iFischer &: Valentil (|2005l ) gives the measured properties of 124 RV planets, 
including those orbiting 45 of the 47 stars listed as planet-bearing in Table 1. The stars left out 
are HD 18445 and and HD 225261. We cannot find any record of these stars having planets, and 
therefore as far as we can tell they are typos in Table 1. 

Since all planets wi th masses above 1 Mjup and periods less than 4 years orbiting stars in the 
Fischer Sz Valenti (j2005l ) list of 850 may be relied upon to have been discovered, we may pick any 
sub-intervals in this range of mass and period, and divide the number of planets falling into these 
intervals by 850 to obtain our normalization. We selected the range 1 -13 M,]„p in mass, and 0. 3-2.5 
AU in semimajor axis. Twenty-eight stars, or 3.29% of the 850 in the lFischer &: Valentil (|2005l ) list, 
have one or more planets in this range. Our inner limit of 0.3 AU excludes the hot J upiters, and thus 



the 3. 29% value provides our final normalization. We note that if we adopt the ICumming et al 
(j2008l ) best-fit power laws, and use the 3.29% normalization to predict the percentage of stars 



having planets with masses between 0.3 and 10 Mj,ip and orbita l perio ds between 2 and 2000 days, 
we find a value of 9.3%, which is close to the ICumming et al.l ()2008l ) value of 10 .5%. The slight 



differe nce is probably not significant, but might be viewed as upward bias in the ICumming et al 



(j2008l ) value due to the inclusion of the hot Jupiters. In any ca se we would not have obtained very 
different constraints if we had used the Cumming et al. (j2008l ) normalization in our Monte Carlo 
simulations. 

For comparison , among the other papers reporting Monte Carlo simulations similar to ours, 
Kasper et al.l ()2007l ) used a normalization of 3% for planets with semimajor axe s of 1-3 AU and 



masse s gre ater than 1 M.i n p. Th is is close to our value of 3.29% for a similar range. iLafreniere et al 



(|2007l ) and iNielsen et al.l (l2008l) fixed a and /3 in their simulations, and l et the normaliza t ion b e 
a fr ee parameter. IChauvin et al.l ()2010l ) obtained their normalization froin ICumming et al.l ()2008l ). 
and Nielsen fc Close ( j2009l ) obtained theirs from Fischer fc Valenti ( 2005 ). 



Juric &: Tremaine I (|2008l ) provide a helpful mathematical description of the eccentricity dis- 



tribution of known RV planets: 



P(6) 



ee 



(5) 
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where P{e) is the probability of a given extrasolar planet's having orbital eccentricity e, e is 
the root of the natural logarithm, and a = 0.3. We find that this mathematical form provides an 
excell ent fit to the distribution of real exoplanet eccentricities from Table 2 of iFischer &: Valenti 
(j2005l ). so we have used it as our probability distribution to generate random eccentricities for the 
Monte Carlo simulations we describe in Section [3] below. 



Constraints on the Distribution of Planets 



3.1. Theoretical Spectra 



Burrows et al.l (j2003l ) present high resolution, flux-calibrated theoretical spectra of giant plan- 



ets or brown dwarfs for ages ranging from 0.1-5.0 Gyr and masses from 1 to 20 Mjup (these are 
available for download from http://www.astro.princeton.edu/~burrows/). We have integrated 
these spectra to give absolute magnitudes in the L' and M filters used in Clio (see Tables [T] and 
12]), and have found that the results can be reasoi iably interpolated to give the L' or M band mag- 
nitudes for all planets of interest for our survey. iBaraffe et al.l (j2003l ) also present models of giant 
planets and brown dwarfs, pre-integrated into magnitudes in the popular infrared bands. These 
models predict slightly better sensitivit y to low mass planets in the L' band and slightly poorer 
sensitivity in the M band, relative to the lBurrows et al.l (120031 ) models. We cannot say if the differ- 
enc e is due to th e sligh tly different filter sets used (MKO for Clio vs. Johnson- Gla ss and Johnson 
for Baraffe et al 



( 2003I')). or if it is intrinsic to the different model spectra used in Burrows et al 



and iBaraffe et al.l (j2003l ). We have chosen to use the lBurrows et al.l (120031) models exclu 



sively herein, to avoid any errors due to the slight filter differences. Since the lBurrows et al.l (j2003l ) 
models predict poorer sensitivity in the L' band, in which the majority of our survey was conducted, 
our decision to use them is conservative. 



3.2. Introducing the Monte Carlo Simulations 



I n common with seve ral other surveys (iKasper et al.l l2007l : iBiller et al.l 120071 : iLafreniere et al 



20071 : IChauvin et al.ll2010l ) we have used our survey null result to set upper limits on planet popula- 
tions via Monte Carlo simulations. In these simulations, we input our sensitivity data in the form of 
tabular files giving the sensitivity in apparent magnitudes as a function of separation in arcseconds 
for each star. Various features of our images could cause the sensitivity at a given separation to vary 
somewhat with position angle: to quantify this, our tabular files give ten different values at each 
separation, corresponding to ten different percentiles ranging from t he worst to th e best sensitivity 
attained at that separation. These files are described in detail in iHeinze et al.l (j2010l ). and are 



available for download from |http : //www . hopewriter . com/Astronomyf iles/Data/SurveyPaper/ 



The Monte Carlo simulations described below allow us to use the observed sensitivity to plan- 
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Table 1. L' Band Absolute Mags from iBurrows et alj (|2003l ) 



Planet Mass 
in Mjup 


Mag at 
0.10 Gyr 


Mag at 
0.32 Gyr 


Mag at 
1.0 Gyr 


Mag at 
3.2 Gyr 


Mag at 
5.0 Gyr 


1.0 


19.074 


23.010 


27.870 


33.50^ 


35.50^ 


2.0 


16.793 


19.351 


23.737 


28.398 


29.479 


5.0 


14.500 


16.397 


18.588 


22.437 


24.407 


7.0 


13.727 


15.390 


17.336 


20.131 


21.574 


10.0 


12.888 


14.437 


16.246 


18.480 


19.466 


15.0 


12.00^ 


13.61'' 


14.773 


16.816 


17.691 


20.0 


11.30^ 


12.98'' 


14.190 


15.967 


16.766 



No models for these very faint planets appear in iBurrows etHI (|2003l '). We have inserted ad hoc values to smooth 
the interpolations. Any effect of the interpolated magnitudes for planets we could actually detect is negligible. 

No models for these bright, hot p lanets appear in burrows et al] l|2003l ). which focuses on cooler objects. We have 
added values from Baraffe et al.1 and then adjusted them to slightly fainter values to ensure smooth interpolations. 



Table 2. M Band Absolute Mags from burrows et al.l \2mi 



Planet Mass 
in Mjup 


Mag at 
0.10 Gyr 


Mag at 
0.32 Gyr 


Mag at 
1.0 Gyr 


Mag at 
3.2 Gyr 


Mag at 
5.0 Gyr 


1.0 


14.974 


16.995 


19.987 


25.0^ 


26.0^ 


2.0 


14.023 


15.313 


17.807 


21.295 


22.163 


5.0 


13.014 


14.017 


15.153 


17.167 


18.537 


7.0 


12.618 


13.561 


14.558 


16.126 


16.909 


10.0 


12.189 


13.096 


14.093 


15.315 


15.951 


15.0 


11.55'' 


12.60'' 


13.370 


14.512 


14.990 


20.0 


11.29'' 


12.21'' 


13.069 


14.122 


14.580 



^No models for these very faint planets appear in IBurrows et al.1 (|2003f ). We have inserted ad hoc values to smooth 
the interpolations. Any effect of the interpolated magnitudes for planets we could actually detect is negligible. 



No models for these bright, hot p l anets appear in IBurrows et all (|2003l ) . which focuses on cooler objects. We 
have added values from Baraffe et al. and then adjusted them to slightly fainter values to ensure smooth 

interpolations. 
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ets in our survey to calculate directly the probability of a given parameter or set of parameters 
describing the exoplanet population. This, in turn, allows us to constrain these parameters at a 
given confidence level. This is a maximum likelihood technique that allows us to incorporate all 
the individual probability functions of the data, as well as parameterized models of the exoplanet 
population. The approach is similar to a Bayesian approach. However, we also use the results of 
the simulations to set confidence limits to the parameters, a more classical approach. 

Each Monte Carlo simulation runs with given planet distribution power law slopes a and /?, and 
a given outer truncation value Rtrunc for the semimajor axis distribution. Using the normalization 
described in Section [21 the probability Ppian of any given star having a planet between 1 and 20 
Mjup is then calculated from the input a, (3, and Rtrunc- In each realization of our survey, each 
star is randomly assigned a number of planets, based on Poisson statistics with mean Ppian- In 
most cases Ppian « 1, so the most likely number of planets is zero. If the star turns out to have 
one or more planets, the mass and semimajor axis of each are random ly selected from the inpu t 



power law distributions. The eccentricity is randomly selected from the lJuric &: Tremaine I (j2008l ) 
distribution, and an inclination is randomly selected from the distribution P{i) gc sm{i). If the star 
is a binary, the planet may be dropped from the simulation at this point if the orbit seems likely 
to be unstable. In general, we consider circumstellar planets to be stable as long as their apastron 
distance is less than 1/3 the projected distance to the companion star, and circumbinary planets 
to be stable as long as their periastron distance is at least three times greater than the projected 
separation of the binary. For planets orbiting low-mass secondaries, a smaller limit on the apastron 
distance is sometimes imposed, while often circumbinary planets required such distant orbits that 
they were simply not considered; the details are given in Table [H For each planet passing the 
orbital stability checkpoint, a full orbit is calculated using a binary star code written by one of 
us (M. K.). The projected separation in arcseconds is found, and the ma gnitude of the planet is 



calculated from its mass, distance, and age using the iBurrows et al.l (120031 ) models. 



Two further random choices complete the determination of whether the simulated planet is 
detected. First, one of the ten percentiles given in the sensitivity files is randomly selected. Com- 
bined with the separation in arcseconds, this selection specifies the sensitivity of our observation 
at the location of the simulated planet. The second random choice is needed because planets ap- 
pearing at low significance in our images would have a less than 100% chance of being confidently 
detected. Our blind sensitivity tests using fake planets placed in our raw data showed that we could 
confirm 97% of lOcr sources, 46% of 7a sources, and 16% of 5a sources, w here cr is a measure of the 



PSF-scale noise in a given region of the image (see iHeinze et al.l (|2010l ) for details). This second 
and final random choice in our Monte Carlo simulations is therefore arranged to ensure that a ran- 
domly selected 16% of planets with 5-7a significance, and 46% of planets with 7-lOcr significance, 
are recorded in the simulation as detected objects. Although we have 97% completeness at lOcr, 
we choose to consider 100% of simulated planets with lOci or greater significance to be detected, 
because at only slightly above lOcr the true completeness certainly becomes 100% for all practical 
purposes. Note that we have conservatively allowed the detection probabilities to increase stepwise. 
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rather than m a continuous curve, from 5 to 10a: that is, in our Monte Carlo simulations, planets 
with 5-7(7 significance are detected at the 5a rate from our blind sensitivity tests, while those with 
7-10(7 significance are detected at the 7a rate. 

The low completeness (16%) at 5a, as determined from our blind sensitivity tests using fake 
planets, may seem surprising. In these tests we distinguished between planets that were suggested 
by a concentration of unusually bright pixels ('Noticed'), or else confidently identified as real sources 
('Confirmed'). Many more planets were noticed than were confirmed: for noticed planets, the 
rates are 100% at 10a, 86% at 7a, and 56% at 5a. However, very many false positives were also 
noticed, so sources that are merely noticed but not confirmed do not represent usable detections. 
The completeness levels we used in our Monte Carlo simulations (16% at 5a and 46% at 7a) 
refer to confirmed sources. No false positives were confirmed in any of our blind tests. Followup 
observations of suspected sources are costly in terms of telescope time, so a detection strategy with 
a low false-positive rate is important. 

Though sensitivity estimators (and therefore the exact mean ing of 5a) differ am ong planet 



imaging surveys, ours was quite conservative, as is explained in iHeinze et al.l (j20ld ). The low 
completeness we find at 5a, which has often been taken as a high-completeness sensitivity limit, 
should serve as a warning to future workers in this field, and an encouragement to establish a 
definitive significance-completeness relation through blind sensitivity tests as we have done. 



Note that our blind sensitivity tests, covered in iHeinze et al.l (|20ld ). are completely distinct 
from the Monte Carlo simulations covered herein. The blind tests involved inserting a little over a 
hundred fake planets into our raw image data to establish our point-source sensitivity. In our Monte 
Carlo work we simulated the orbits, masses, and brightnesses of millions of planets, and compared 
them to our previously-established sensitivity limits to see which planets our survey could have 
detected. 



3.3. A Detailed Look at a Monte Carlo Simulation 

To evaluate the significance of our survey and provide some gu idance for future work, we have 



analyzed in detail a single Monte Carlo simulation. We chose the ICumming et al.l (|2008l ) best fit 
values of (3! = —1.31 and (3 = —0.61, with the semimajor axis truncation radius set to 100 AU. 
Planets could range in mass from 1 to 20 Mjup. As described in Section [2] above, we normalized the 
planet distributions so that each star had a 3.29% probability of having a planet with semimajor 
axis between 0.3 and 2.5 AU and mass between 1 and 13 Mj^p. The simulation consisted of 50,000 
realizations of our survey with these parameters. In all, 505,884 planets were simulated, of which 
51,879 were detected. 

In 38% of the 50,000 realizations, our survey found zero planets, while 37% of the time it 
found one, and 25% of the time it found two or more. The planet distribution we considered in 
this simulation cannot be ruled out by our survey, since a null result such as we actually obtained 
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turns out not to be very improbable. 

The large number of survey realizations in our simulation allows the calculation of precise 
statistics for potentially detectable planets. The median mass of detected planets in our simulation 
was 11.36 Mjup, the median semimajor axis was 43.5 AU, the median angular separation was 2.86 
arcsec, and the median significance was 21Aa. This last number is interesting because it suggests 
that, for our survey, any real planet detected was likely to appear at high significance, obvious 
even on a preliminary, 'quick-look' reduction of the data. This suggests that performing such 
reductions at the telescope should be a high priority, to allow immediate confirmation and followup 
if a candidate is seen. Figure [1] presents as a histogram the significance of all planets detected in 
this Monte Carlo simulation. 

We suspected that there would be a detection bias toward very eccentric planets, because these 
would spend most of their orbits near apastron, where they would be easier to detect. This bias 
did not appear at any measurable level in our simulation. However, there was a weak but clear 
bias toward planets in low-inclination orbits, which, of course, spend more of their time at large 
separations from their stars than do planets with nearly edge-on orbits. 

A concern with any planet imaging survey is how strongly the results hinge on the best (i.e. 
nearest and youngest) few stars. A survey of 54 stars may have far less statistical power than the 
number would imply if the best two or three stars had most of the probabilty of hosting detectable 
planets. Table [3] gives the percentage of planets detected around each star in our sample based on 
our detailed Monte Carlo simulation. Due to poor data quality, binary orbit constraints, or other 
issues, a few stars had zero probability of detected planets given the distribution used here. In 
general, however, the likelihood of hosting detectable planets is fairly well distributed. 

In Table m we give the details of planetary orbital constraints used in our Monte Carlo simula- 
tions for each binary star we observed, complete with the separations we measured for the binaries. 
Note that HD 96064 B is a close binary star in its own right, so planets orbiting it were limited 
in two ways: the apastron could not be too far out, or the orbit would be rendered unstable by 
proximity to HD 96064 A - but the periastron also could not be too far in, or the binary orbit of 
HD 96064 Ba and HD 96064 Bb would render it unstable. Planets individually orbiting HD 96064 
Ba or HD 96064 Bb were not considered in our survey, since to be stable the planets would have 
to be far too close-in for us to detect them. The constraints described in Table |4] account for most 
of the stars in Table [3] with few or no detections reported. 

A final question our detailed simulation can address is how important the M band observations 
were to the survey results. In Tabled we show that when M band observations were made, they 
did substantially increase the number of simulated planets detected. 
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Fig. 1. — Histogram of detection si gnificance for th e 51,87 9 simulated planets detected in 50,000 
realizations of our survey with the I Gumming et alJ (120081 ) distribution (a = —1.31, /? = —0.61) 
truncated at 100 AU. Our detection rates went down for significance less than 10 a, but some 5-7 
cr planets are still detected. The relatively high median significance of 21.4 a suggests any detected 
planet would most likely be quite obvious - a good argument for doing 'quick-look' data reductions 
as soon as possible at the telescope. 
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Table 3. Percentage of Detected Planets Found Around Each Star 





% of Total 


A/Tpr] li^n 
X V i vj. ± cx J. J. 


IVTon 1 f\ Ti 

J. V J. Vy J. J. (. Ij i i. 


A Tpr] 1 Pi Tl 






Mass 




'-J 1-' C*iX (Ai \J i\Jl.± 


GJ 117 


6.07 


7.66 Mjup 


39.36 AU 


3.64 arcsec 


e Eri 


5.83 


6.98 Mjup 


18.26 AU 


4.35 arcsec 


HD 29391 


5.80 


8.14 Mjup 


49.13 AU 


2.71 arcsec 


GJ 519 


4.74 


10.44 Mjup 


40.51 AU 


3.28 arcsec 


GJ 625 


4.67 


9.72 Mjup 


29.18 AU 


3.48 arcsec 


GJ 5 


4.45 


9.60 Mjup 


53.42 AU 


3.08 arcsec 


BD+60 1417 


3.95 


11.58 Mjup 


44.48 AU 


2.05 arcsec 


GJ 355 


3.81 


9.71 Mjup 


53.91 AU 


2.34 arcsec 


GJ 354.1 A 


3.67 


9.58 Mjup 


60.12 AU 


2.64 arcsec 


GJ 159 


3.57 


9.73 Mj^p 


57.95 AU 


2.71 arcsec 


GJ 349 


3.35 


11.38 Mjup 


44.40 AU 


3.17 arcsec 


61 Cyg B 


3.29 


11.32 Mjup 


19.53 AU 


4.08 arcsec 


GJ 879 


3.03 


11.18 Mjup 


36.84 AU 


3.69 arcsec 


GJ 564 


2.94 


10.67 Mjup 


56.80 AU 


2.70 arcsec 


GJ 410 


2.93 


12.78 Mjup 


41.83 AU 


3.03 arcsec 


GJ 450 


2.89 


12.90 Mjup 


38.72 AU 


3.66 arcsec 


GJ 3860 


2.68 


12.70 Mjup 


49.72 AU 


2.69 arcsec 


HD 78141 


2.58 


12.47 Mjup 


57.00 AU 


2.24 arcsec 


BD+20 1790 


2.51 


12.14 Mjup 


58.33 AU 


2.02 arcsec 


GJ 278 C 


2.20 


12.68 Mjup 


54.56 AU 


3.04 arcsec 


GJ 311 


2.19 


12.55 Mjup 


52.07 AU 


3.20 arcsec 


HD 113449 


2.17 


12.52 Mj,p 


59.31 AU 


2.29 arcsec 


GJ 211 


2.10 


13.59 Mjup 


50.51 AU 


3.30 arcsec 


BD+48 3686 


2.08 


12.56 Mjup 


55.05 AU 


2.01 arcsec 


GJ 282 A 


2.05 


13.39 Mjup 


49.85 AU 


2.99 arcsec 


GJ 216 A 


2.03 


12.71 Mjup 


42.98 AU 


4.21 arcsec 


61 Cyg A 


1.97 


13.70 Mjup 


20.94 AU 


4.54 arcsec 


HD 1405 


1.54 


13.13 Mjup 


66.34 AU 


2.04 arcsec 


HD 220140 A 


1.54 


11.73 Mj„p 


36.85 AU 


1.73 arcsec 


HD 96064 A 


1.49 


12.63 Mjup 


46.64 AU 


1.75 arcsec 


HD 139813 


1.43 


14.33 Mjup 


59.71 AU 


2.37 arcsec 


GJ 380 


0.92 


15.76 Mjup 


25.31 AU 


4.21 arcsec 
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3.4. Monte Carlo Simulations: Constraining the Power Laws 

The planet distribution we used in tlie single Monte Carlo simulation described above could 
not be ruled out by our survey. To find out what distributions could be ruled out, we performed 
Monte Carlo simulations assuming a large number of different possible distributions, parametrized 
by the two power law slopes a and /3, and by the outer semimajor axis truncation radius Rtrunc- 
R egardless of the va l ues of a and /?, each simulation was normalized to match the RV statistics 



of iFischer Valentil (|2005l ): any given star had 3.29% probability of hosting a planet with mass 
between 1 and 13 Mj^p and semimajor axis between 0.3 and 2.5 AU. The mass range for simulated 
planets was 1-20 Mj^p. 

We tested three different values of a: -1.1, -1.31, and -1.51, roughly correspon ding to the most 



optini istic permitted, the best fit, and the most pessimistic permitted values from lCumming et al 



(|2008l ). For each value of a, we ran simulations spanning a wide grid in terms of /3 and Rtmnc- In 
constrast to the extensive results described in Section [3.31 the only data saved for these simulations 
was the probability of finding zero planets. Since we did in fact obtain a null result, distributions 
for which the probability of this was sufficiently low can be ruled out. 

Figures [2] and [3] show the probability of a null result as a function of /3 and Rtmnc for our 
three different values of a. Figure [2] presents constraints based on a = —1.31, the best-fit value 
from RV statistics, while Figure [3] compares the optimistic case a = —1.1 and the pessimistic 
case a = —1.51. Each pixel in these figures represents a Monte Carlo simulation involving 15,000 
realizations of our survey; generating the figures took several tens of hours on a fast PC. Contours 
are overlaid at selected probability levels. Regions within the 1%, 5%, and 10% contours can, of 
course, be ruled out at the 99%, 95%, and 90% con fidence levels respective ly. For example, we find 



that the most optimistic power laws allowed by the lCumming et al.l (j2008l ) RV statistics, a = —1.1 
and /3 = —0.46, are ruled out with 90% confidence if Rtmnc is 110 AU or greater. Similarly, 
a = —1.51 and j3 = —0.3, truncated at 100 AU, is ruled out. Though /3 = 0.0 is not physically 
plausible, previous work has sometimes used it as an example: for a = —1.31, we rule out j3 = 0.0 
unless Rtmnc is less than 38 AU. 



3.5. Model-independent Constraints 



It is also possible to place constraints on the distribution of planets without assuming a power 
law or any other particular model for the statistics of planetary masses and orbits. Note well that 
by "model-independent" in this context, we mean independent only of models for the statistical 
distributions of planet s in terms of M and a - not independent of models of planetary spectra such 
those we obtain from iBurrows et al.l (|2003l ). The latter are our only means of converting from 
planetary mass and age to detectable flux, and as such they remain indispensable. 



To place our model-independent constraints, we performed an additional series of Monte Carlo 
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Table 3 — Continued 





% of Total 


Median 


Median 


Median 


Star Name 


Detected Planets 


Mass 


Semimaior Axis 


Separation 


GJ 896 A 


0.61 


12.43 Mjup 


6.47 AU 


0.98 arcsec 


GJ 860 A 


0.38 


11.58 Mjup 


53.26 AU 


6.62 arcsec 


r Ceti 


0.38 


17.19 Mjup 


25.49 AU 


5.52 arcsec 


GJ 896 B 


0.34 


11.40 Mj,p 


6.78 AU 


1.14 arcsec 


^ Boo B 


0.32 


12.07 Mjup 


8.25 AU 


1.36 arcsec 


HD 220140 B 


0.28 


12.04 Mjup 


25.92 AU 


1.37 arcsec 


C Boo A 


0.24 


12.89 Mjup 


8.72 AU 


1.50 arcsec 


GJ 659 B 


0.21 


17.71 Mjup 


62.54 AU 


2.81 arcsec 


GJ 166 B 


0.17 


16.12 Mjup 


6.19 AU 


1.34 arcsec 


GJ 684 A 


0.17 


14.93 Mjup 


85.98 AU 


4.87 arcsec 


HD 96064 B 


0.13 


14.43 Mjup 


38.55 AU 


1.60 arcsec 


GJ 505 B 


0.12 


15.94 Mjup 


17.11 AU 


1.61 arcsec 


GJ 166 C 


0.10 


15.56 Mjup 


6.43 AU 


1.52 arcsec 


GJ 505 A 


0.07 


16.32 Mjup 


18.08 AU 


1.75 arcsec 


GJ 702 A 


0.02 


15.90 Mjup 


6.21 AU 


1.50 arcsec 


GJ 684 B 


None 


NA 


NA 


NA 


GJ 860 B 


None 


NA 


NA 


NA 


GJ 702 B 


None 


NA 


NA 


NA 


HD 77407 A 


None 


NA 


NA 


NA 


GJ 659 A 


None 


NA 


NA 


NA 


GJ 3876 


None 


NA 


NA 


NA 


HD 77407 B 


None 


NA 


NA 


NA 



Note. — This table applies to our detailed Monte Carlo simulation with 50,000 
survey realizations run using a = —1.31, /3 = —0.61, and semimajor axis trunca- 
tion radius 100 AU. Of all the simulated planets that were detected, we present 
here the percentage that were found around each given star, and the median mass, 
semimajor axis, and projected separation for simulated planets found around each 
star. The table thus indicates around which stars our survey had the highest like- 
lihood of detecting a planet. Many stars with poor likelihood are binaries, with 
few stable planetary orbits possible. 
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Table 4. Constraints on Simulated Planet Orbits Around Binary Stars 







constraints on 


constraints on 


constraints on 




separation 


circumprimary 


circumsecondary 


circumbinary 


Star Name 


(arcsec) 


apastron 


apastron 


periastron 


HD 220140 AB 


10.828 


<3.61 asec (71.3 AU) 


<2.17 asec (42.8 AU) 


No Stable Orbits 


HD 96064 AB 


11.628 


<3.88 asec (95.6 AU) 


<2.33 asec (57.3 AU) 


No Stable Orbits 


HD 96064 Bab 


0.217 


No Stable Orbits 


No Stable Orbits 


>0.65 asec (16.1 AU) 


GJ 896 AB 


5.366 


<1.79 asec (11.8 AU) 


<1.79 asec (11.8 AU) 


No Stable Orbits 


GJ 860 AB 


2.386 


<0.79 asec (3.17 AU) 


<0.60 asec (2.41 AU) 


>7.15 asec (28.7 AU) 


^ Boo AB 


6.345 


<2.12 asec (14.2 AU) 


<2.12 asec (14.2 AU) 


No Stable Orbits 


GJ 166 BC 


8.781 


<2.20 asec (10.6 AU) 


<2.20 asec (10.6 AU) 


No Stable Orbits 


GJ 684 AB 


1.344 


<0.45 asec (6.34 AU) 


<0.27 asec (3.80 AU) 


>4.03 asec (56.8 AU) 


GJ 505 AB 


7.512 


<2.50 asec (29.8 AU) 


<2.50 asec (29.8 AU) 


No Stable Orbits 


GJ 702 A 


5.160 


<1.76 asec (8.85 AU) 


<1.32 asec (6.64 AU) 


>15.9 asec (79.7 AU) 


HD 77407 AB 


1.698 


<0.57 asec (17.2 AU) 


<0.34 asec (10.2 AU) 


>5.11 asec (153.7 AU) 



Note. — Planets orbiting the primary in a binary star were considered to be de-stabilized by the 

gravity of the secondary if their apastron distance from the primary was too large. Similarly, planets 
orbiting the secondary had to have small enough apastron distances to avoid being de-stabilized by 
the primary. Circumbinary planets had to have a large enough periastron distance to avoid being de- 
stabilized by the differing gravitation of the two components of the binary. Note that HD 96064B is itself 
a tight binary star, so planets orbiting it had both a minimum periastron and a maximum apastron. 
Constraints are given in AU as well as arcseconds so constraints can easily be compared with actual or 
hypothetical planetary systems. 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 
Semimajor Axis Power Law Slope 



Fig. 2. — Probability of our survey detecting zero planets, as a function of the power law slope of the 
semimajor axis distribution /?, where ^ oc , and the outer truncation radius of the semimajor axis 
distribution. Here, the slope of the mass distribution a has been taken as -1.31, where oc M". 
Since we found no planets, distributions that lead to a probability P of finding no planets are ruled 
out at the 1 — P confidence level: for example, the region above and to the right of the 0.1 contour 
is ruled out at the 90% confidence level 
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Table 5. Importance of the M Band Data 





Total simulated 


2-band 


L'-only 


M-only 


Star Name 


detections 


detections 


detections 


detections 


e Eri 


2850 


46.98% 


8.28% 


44.74% 


61 Cyg B 


1610 


52.73% 


1.55% 


45.71% 


61 Cyg A 


965 


63.01% 


22.80% 


14.20% 


^ Boo B 


157 


61.15% 


18.47% 


20.38% 


^ Boo A 


115 


60.00% 


18.26% 


21.74% 


GJ 702 A 


9 


22.22% 


0.00% 


77.78% 



Note. — The usefulness of M band observations, based on our 
detailed Monte Carlo simulation. When M band observations were 
made of a given star, they did substantially increase the number of 
simulated planets detected around that star. 




-1.0 -0.8 -0.6 -0.4 -0.2 0.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 

Semimajor Axis Power Law Slope Semimajor Axis Power Law Slope 



Fig. 3. — Probability of our survey detecting zero planets, as a function of the power law slope of 
the semimajor axis distribution /?, where ^ oc a^, and the outer truncation radius of the semimajor 
axis distribution. Here, the slope of the mass distribution a has been taken as -1.1 (left) and -1.51 
(right), where tx M". Since we found no planets, distributions that lead to a probability P of 
finding no planets are ruled out at the 1 — P confidence level: for example, the regions above and 
to the right of the 0.1 contours are ruled out at the 90% confidence level 
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simulations on a grid of planet mass and orbital semimajor axis. For each grid point we seek 
to determine a number P{M,a) such that, with some specified level of confidence (e.g., 90%), 
the probability of a star like those in our sample having a planet with the specified mass M and 
semimajor axis a is no more than P{M,a). We determine P{M,a) by a search: first a guess is 
made, and a Monte Carlo simulation assuming this probability is performed. If more than 10% of 
the realizations of our survey turn up a null result, the guessed probability is too low; if less than 
10% turn up a null result, the probability is too high. It is adjusted in steps of ever-decreasing size 
until the correct value is reached. 

Figure H] shows the 90% confidence upper limit on P{M, a) as a function of mass M and 
semimajor axis a. Each pixel represents thousands of realizations of our survey, with P{M, a) 
finely adjusted to reach the correct value. Contours are overplotted showing where P{M, a) is 
less than 8%, 10%, 25%, 50%, and 75%, with 90% confidence. Note that P{M,a), the value 
constrained by our simulations, is a probability rather than a fixed fraction. The probability is 
the more scientifically interesting number, but is harder to constrain. For example, if 3.7% is the 
fraction of the actual stars in our sample that have planets with easy-to-detect properties, there 
are 2 such planets represented in our 54-star survey. However, if the probability of a star like those 
in our sample having such a planet is 3.7%, there is still a nonzero probablity (13% in this case) 
that no star in our sample actually has such a planet. 

The results presented in Figure|4]can be interpreted as model-independent constraints on planet 
populations. For example, with 90% confidence we find that less than 50% of stars with properties 
like those in our survey have a 5 Mjup or more massive planet in an orbit with a semimajor axis 
between 30 and 94 AU. Less than 25% of stars like those in our survey have a 7 Mjup or more 
massive planet between 25 and 100 AU, less than 15% have a 10 Mjup or more massive planet 
between 22 and 100 AU, and less than 12% have a 15 Mjup or more massive planet/brown dwarf 
between 15 and 100 AU. Going to the most massive objects considered in our simulations, we can 
set limits ranging inward past 10 AU: we find that less than 25% of stars like those surveyed have 
a 20 Mjup object orbiting between 8 and 100 AU. These constraints hold independently of how 
planets are distributed in terms of their masses and semimajor axes. 

HR 8799 appears to have a remarkable system of three massive planets, seen at projected 
dista nces of 24, 38, and 68 AU, with masses of roughly 10, 10, and 7 Mjup, respectively (jMarois et al 



20081 ). Using a Monte Carlo simulation like those used to create Figure [U we find with 90% 
confidence that less than 8.1% of stars like those in our survey have a clone of the HR 8799 planetary 
system. For purposes of this simulation we adopted the masses above, and set the planets' orbital 
radii equal to their projected separations. Our 8.1% limit represents a step toward determining 
whether or not systems of massive planets in wide orbits are more common around more massive 
stars such as HR 8799 than FGK stars such as those we have surveyed. 
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Fig. 4. — 90% confidence level upper limits on the probability P{M, a) that a star like those in our 
survey will have a planet of mass M and semimajor axis a. This plot shows, for example, that our 
survey constrains the abundance of 10 Mjup or more massive planets with orbital semimajor axes 
between 22 and 100 AU to be less than 15% around sun-like stars. The abundance of 5 Mj^p or 
more massive planets between 25 and 94 AU is constrained to be less than 50%. The latter range 
does not extend all the way to 100 AU because our sensitivity to planets in very distant orbits 
decreases somewhat due to the possibility of their lying beyond our field of view. 
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3.6. Our Survey in the Big Picture 



The surveys of iKasper et al.l (j2007l ) and iBiller et alj (j2007l ) , have set c onstraints on the dis - 
tributions of e xtrasolar planets simil ar to those we present herein, while iNielsen et al.l (j2008l ) 
and especiallv iLafreniere et al.1 (1200711 hav e set stronger constraints. More recent analyses by 



Nielsen &: Closd (120091 ) and IChauvin et alJ (l20ld l also provide constraints on th e planetary dis- 



tribut ion. For example. iNielsen &: Closd (120091 ) provide a 68% confidence that the lCumming et al 
(j2008l ) distribution can be excluded for a truncation radius of 28 AU. However, if different models 
are used this number jumps to 83 AU. Chauvin et al. indicate a similar limit from analyzing their 
results using Baraffe et al. 2003 models. For the standard parameters they indicate a maximum 
permitted truncation radius of approximately 35 AU. In this context, the results presented here 
provide looser constraints on the planet distribution, but provide an independent check on the 
model-dependent systematic errors which may exist with shorter wavelength data, due to incorrect 
model brightness estimates or age determination. 

Theoretical spectra of self-luminous extrasolar planets are very poorly constrained observation- 
ally. The recent d etections of possible planets aroun d HR 8799 (IMarois et al.l 120081 ). Fomalhaut 
(jKalas et al.ll2008l ). and /3 Pic (jLagrange et al.l l2009l ) are either single-band (/3 Pic) or only be- 
ginning to be evaluated at multiple wavelengths (HR 8799, Fomalhaut). The candidate planets 
orbiting HR 8799 and (3 Pic are hotter than we would expect to find orbiting middle-aged stars 
such as those in our survey, while HST photometry of Fomalhaut b suggests much of its brightness 
is starlight reflected from a circumplanetary dust disk. Our survey, and other exoplanet surveys, 
must therefore be interpreted using models of planetary spectra that are not yet well-tested against 
observations. 

Such models predict brightnesses in the H band, and particularly in narrow spectral win- 
dows within the H band, that are enormously in excess of black body fluxes. The constraints 



set by 


Masciadri et al. 


Nielsen & Close ( 


2009) 



of remarkable brightness in the H band. The V and M bands that we have used are nearer 
the blackbody peaks of low-temperature self-luminous planets, and might be expected to be more 
reliable. 



However, iLeggett et al.l (|2007j) and iReid &: Crua (|2002l ) suggest that the M b a nd br ightness 
at least of hotter extrasolar planets will be less than predicte d by Burrows et al. ( 20031 ) due to 
above-equilibrium concentrations of CO from convective mixing. iHubenv &: Burrowsl (j2007l ) present 
new models indicating the effect is present for planets with Tefj ranging from 600 to 1800K. The 
maximum M band flux su pression is about 40% , and flux supression disappears completely for 



Teff below 500K. Based on iBurrows et al.l (j2003l ). this Teg value corresponds to planets of about 
3.5, 6.5, 12, and 15 Mjup at ages of 100 Myr, 300 Myr, 1 Gyr, and 2 Gyr, respectively. In many 
cases our M band observations were sensitive to planets at lower masses than these values, and 
therefore Teff lower than 500K, implying that the CO flux supression would have no effect on our 
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mass limits. In other cases our M band sensitivity did not extend so low. However, given that 
M band observations formed a relatively small part of our survey, and CO supression would affect 
only a fraction even of them, the total effect on the statistical conclusions of our survey should be 
entirely negligible. 



Theoretical spectra such as those of iBurrows et al.l ([2003) may or may not be more reliable in 
the L' and M bands than at shorter wavelengths. However, so long as the models remain poorly 
constrained by observations at every wavelength, conclusio ns based on observat ions at multiple 
wavelengths will be more secure. Our survey, with that of iKasper et al.l ( 20071 ). has diversified 
planet imaging surveys across a broader range of wavelengths. 



In another sense our survey differs even from that of lKasper et al.l (120071 ) : we have investigated 
older stars. This is significant because planetary systems up to ages of several hu ndred Myr may still 



be undergoing substantia l dynamical evolution due to planet-planet interactions (jJuric &: Tremaine 



20081 : iGomes et al.ll2005l'). Our survey did not necessarily probe t he same planet population as, for 
example, those of lKasper et al.l (|2007l ) and IChauvin et al.l (j2010l ). 



Finally, theoretical models of older planets are likely more reliable than for younger ones 
as these planets are further from their unknown starting conditions and moving t oward a well- 



under stood, stable configuratiori such as that of Jupiter. It has been suggested b y iMarlev et al 



200 



3) and 



Fortnev et a 



(120081 ) that theoretical planet models such as those of iBurrows et al 



2OO3I ) and tearaffe et all (|2003l ) may overpredict the brightness of young (< 100 Myr) planets by 



orders of magnitude, while for older planets the models are more accurate. 

We have focused on nearby, mature star systems, and have conservatively handled the ages 
of stars. This makes our survey uniquely able to confirm that the rarity of giant planets at large 
separations around solar-type stars, first noticed in surveys strongly weighted toward young stars, 
persists at older system ages. It is not an artifact of model inaccuracy at young ages due to unknown 
initial conditions. 



4. The Future of the L' and M Bands 



In the L' and M bands, the sky bri ghtness is much worse than at shorter wavelengths. How- 
ever, models (e.g., IBurrows et al.l (j2003l )) predict that in the L' and M bands, planets fade less 
severely with increasing age (or, equivalently, decreasing Tes). Also, planet/star flux ratios are 
more favorable in the L' and M bands than at shorter wavelengths such as the H and Ks bands. 

It makes sense to use the L' and M bands o n bright stars , wher e the planet/star flux ratio is a 
more limiting factor than the sky brightness. In lHeinze et al.l (|2008l ). we have shown that M band 
observations tend to do better than those at shorter wavelengths at small separations from bright 
stars. 



The L' and M bands are most useful, however, for detecting the lowest temperature planets. 
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which have the reddest H — L' and H — M colors. Such very low temperature planets can only be 
detected around the nearest stars, so it is for very nearby stars that L' and M band observations 
are most useful. For distant stars, around which only relativly high Teg planets can be detected, 
the H and Ks bands are much better. We will now quantitatively describe the advantage of L' 
and M band observations over shorter wavelengths for planet-search observations of nearby stars. 

Most AO planet searches to date have used the H and Ks bands, or specialized filters in 
the same wavelength regime. While the Ks band has been used extensively to search for planets 
around young stars (jMasciadri et al.l l2005l : IChauvin et al.l l20ld ) , our comparison here will focus 
on the H band regime. Model s indicate it offers better sensitivity tha n the Ks band except for 
planets younger than 100 Myr (jBurrows et al.l l2003l : iBaraffe et al.ll2003l ). and most of the stars we 
will suggest the L' and the M bands are useful for will be older than this. The m ost sensitive H- 
regime planet search observations made to date are those of lLafreniere et al.l (120071 ) , in part because 
of their optimized narrow-band filter. They attained an ef fective background-lim i ted point-source 
sensit ivity of about H = 23.0. Based on the models of iBurrows et al.l (j2003l ). iLafreniere et al 



(|2007l ) would have set better planetary mass limits than our observations around all of our own 
survey targets except the very nearest objects, such as e Eri and 61 Cyg. Thus, at present, the 
iJ-regime delivers far better planet detection prospects than the L' and M bands for most stars. 

However, as detector technology improves, larger telescopes are built, and longer planet de- 
tection exposures are attempted, the sensitivity at all wavelengths will increase. This means that 
low-temperature planets, with their red IR colors, will be detectable at larger distances, and the 
utility of the L' and especially the M bands will increase. In Figure [5] we show the minimum 
detectable planet mass for hypothetical stars at 10 and 25 pc distance as a function of the increase 
over current sensitivity in the H, L', and M bands, and in Figure [6] we present t he same compar- 
i son f or a star at 5 pc. We have taken current sensitivity to be -ff = 23.0 (i.e., ILafreniere et al. 
(120071)'). L' = 16.5, a. nd M = 13.5 (i.e., the present work, scaled to an 8m telescope such as 



Lafreniere et al 



(120071) used). These are background limits, not applicable close to bright stars. 



Based on Heinze et al 



(I2OO8), we believe the L' and M bands will do even better relative to H 
closer to the star where observations are no longer background limited. Of course H band obser- 
vations with next-generation extreme AO systems such as GPI and SPHERE w ill offer improved 



perfor mance close to the star, but advances in M-band AO coronography (e.g. iKenworthy et al 



(|2007l )). will also improve the longer- wavelength results. In any case. Figures [5] and [6] compare 
background-limited performance only. 



The supressio n of flux in the M band due to elevated levels of CO (jLeggett et al 



2007 



Reid &: Cruz 120021') does no t appl y to planets at the low temperatures relevant for Figures [5] and 



El Based on I Bur rows et al.l (120031 ) . the entire mass range covered by both Figures corresponds to 
planets with Tgfj below 500K, except for planets with masses above 6.5 Mjup in the left panel of 
Figure [5] (25 pc distance, 300 Myr age). This upper section of the 25 pc, 300 My r pane l is irrelevant 
to the important implications of the figure. According to iHubeny &: Burrowsl (120071 ). there is no 
supression of the M band for effective temperatures below 500K. 
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Fig. 5. — Minimum detectable planet mass in units of Mjup for stars at 25pc (left) and lOpc (right), 
in the H, L' , and M bands, as a function of increase over current sensitivity. We have taken current 
sensitivities to he H = 23.0, L' = 16.5, and M = 13.5. While the H band will likely remain the 
wavelength of choice for planet search observations of stars at 25 pc and beyond, an increase of 
only 2.4 mag over current sensitivities, even though paralleled by an equal increase in H band 
sensitivity, will render the M band more sensitive than H for planets around all stars nearer than 
10 pc. The relative effectiveness of different wavelengths depends sensitively on the distance to a 
star system, but it is essentially independent of the stellar age, as explained in the text. 
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Fig. 6. — Minimum detectable planet mass in units of Mj^p for stars at 5pc, in the H, L' , and M 
bands, as a function of increase over current sensitivity. We have taken current sensitivities to be 
H = 23.0, L' = 16.5, and M = 13.5. Given only a 1 magnitude increase in M band sensitivity, 
paralleled by an equal increase at H band, the M band would be the best wavelength for planet 
search observations around all stars nearer than 5 pc. While the sensitivity increases required 
to render M preferable in Figure [5] require substantial improvements to existing instruments and 
telescopes, the 1 mag increase required at 5 pc could be obtained by simply increasing the exposure 
time. As with Figure [5l this result concerning the relative effectiveness of different wavelengths is 
independent of stellar age, to first order. 
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We have deliberately chosen the characteristics of the hypothetical stars in Figures [5] and [6] 
to be less good than the best available planet search candidates, so that in each case stars closer 
and/or younger than the example actually exist. Using the very youngest stars would also have 



resulted in sensitivities better than 1 Mjup , a mass regime not covered by the I Burrows et alj (|2003l ) 
models used in the Figures. 

Figures [5] and [6] illustrate three very important points. First, the L' band appears to have only 
secondary usefulness since either the H band or the M band always offers sensitivity to lower-mass 
planets. Second, Figure [6] shows that with a relatively minor increase of 1 magnitude in sensitivity, 
the M band will be sensitive to lower-mass planets around all stars within 5 pc than can be detected 
with H band observations, even if the H band sensitivity increases the same amount. Third, Figure 
[5] shows that the advantage of the M band decreases with increasing distance, but that as larger 
telescopes and longer exposures increase sensitivities to 2.5 mag above present levels, the M band 
will be superior to H out to 10 pc. With an increase of 4 mag, the M band would surpass H 
out to 25 pc - but as such a large sensitivity increase would be difficult to achieve, H band will 
likely remain the primary wavelength for stars at 25 pc and beyond. For stars closer than 10 pc, 
however, the M band already offers excellent sensitivity that has barely been exploited so far. 
Given reasonable sensitivity increases, M should become the primary band for planet searches 
around stars at a distance of 10 pc or less. 

Interesting ly, the c onclu sions of Figures [5] and [6] are essentially independent of age: extensive 
calculations by Heinze (j2007l ) showed that the relative usefulness of different wavelengths had only 



a weak dependence on age, for stars at a fix ed distance - and even this weak age dependence coul d 



change sign on switching from the models of lBurrows et al.l (|2003l ) to those of lBaraffe et al.l (j2003l ). 



This means that if we change the ages of the stars in Figures [5] and [6] but leave the distances the 
same, the L', M, and H band curves will slide up or down but remain essentially fixed in their 
relative positions. For example, given a 3 magnitude increase in sensitivity at both wavelengths, M 
band observations will detect lower mass planets than ff-band ones around a star at 10 pc, whether 
the stellar age is 5 Gyr, 1 Gyr, or 100 Myr. This is to be expected, since if one dials down the 
age of a given hypothetical star system, the Tgfj (and therefore IR color) of the faintest detectable 
planets will remain about the same, though their masses will decrease. 

Again, Figures [5] and [6] apply only to background- limited sensitivity. However, given the much 
more favorable planet/star flux ratios in the M band relative to we would expect the longer 
wavelength observations to remain equally competitive closer to the star. Advances in M band 
coronography will likely parallel the development of H band extreme AO systems such as GPI and 
SPHERE. Though at present they are surpassed in sensitivity by //-regime observations for all but 
the nearest stars, the L' and especially the M bands hold considerable promise for the future. 
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Conclusion 



We have surveyed unusually nearby, mature star systems for extrasolar planets in the L' and 
M bands using the Clio camera with the MMT AO system. By extensive u se of blind sensitivi ty 
tests involving fake planets inserted into our raw data (reported in detail in lHeinze et alj (|20ld )). 
we established a definitive significance vs. completeness relation for planets in our data, which we 
then used in Monte Carlo simulations to constrain planet distributions. 

We set interesting limits on the masses of planets and brown dwarfs in the star systems we 
surveyed, but we did not detect any planets. Based on this null result, we place constraints on 
the power laws that may describe the distribution of extrasolar planets in mass and semimajor 
axis. We also place constraints on planet abundances independent of t he distributions. If the 
dist ribution of pla.nets is a power law with dN oc M'^a^ dMda, the work of ICumming et al.l ( 20081 ) 
and lButler et al.l (|2006l ) indicates that the most optimistic (i.e. planet-rich) case permitted by the 
statistics of known RV planets correponds to about a = —1.1 and /3 = —0.46. Normalizing the 
distribution to be consistent with RV statistics, we find that these values of a and /3 are ruled out 
at the 90% confidence level, unless the semimajor axis distribution is truncated at a radius Rtrunc 
less than 110 AU. Though (3 = 0.0 is not physically plausible, previou s work has some t imes used 
it an example: for a = —1.31, corresponding to the best-fit value from ICumming et al.l (|2008l ). we 
rule out (3 = 0.0 unless Rtrunc is less than 38 AU. Independent of distribution models, with 90% 
confidence no more than 50% of stars like those in our survey have a 5 Mj^p or more massive planet 
orbiting between 30 and 94 AU, no more than 15% have a 10 Mjup planet orbiting between 22 and 
100 AU, and no more than 25% have a 20 Mjup object orbiting between 8 and 100 AU. 

(I2OO7I) and 



Our constrain ts on planet abundances are similar to those pla ced bv lKasper et al 



Biller et al.l (|2007l ). but less tight than those o f iNielsen et al.l (l2008l) and especi allv lLafreniere et al 
(I2OO3), The recent work of iNielsen &: Closel (12002) and Ichauvin et al.l (|20ld ) also placed tighter 
constraints on exoplanet distributions than our survey. However, we have surveyed a more nearby, 
older set of stars than any previous survey, and have therefore placed constraints on a more mature 
population of planets. Also, we have confirmed that a paucity of giant planets at large separations 
from sun-like stars is robustly observed at a wide range of wavelengths. 



The best current H regime observations, those of iLafreniere et al.l (|2007l ). would attain sensi- 
tivity to lower mass planets than did our L' and M band observations for all of our survey targets 
except those lying within 4 pc of the Sun. However, as larger telescopes are built and longer expo- 
sures are attempted, the sensitivity of M band observations may be expected to increase at least as 
fast as that of H band observations (in part because M band detectors are currently a less mature 
technology). As shown in Figures [5] and El a modest increase from current sensitivity levels, even 
if paralleled by an equal increase in H band sensitivity, would render the M band the wavelength 
of choice for extrasolar planet searches around a large number of nearby stars. 
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